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The Equivalent Random method is used for engineering many of 
the telephone overflow-networks in the Bell System. But since this 
method is not directly applicable to the analysis of graded-multiple 
trunk-groups which carry overflow traffic, we extend the method to 
cover such arrangements. The key to this extension is a technique for 
taking correlation into account when combining dependent streams 
of traffic which are themselves more variable than Poisson, In prin- 
ciple, the technique is applicable wherever a stream of overflow traffic 
is divided, submitted to independent trunk groups, and then recom- 
bined. 

The extended Equivalent Random method provides adequate esti- 
mates of load-service relations for graded multiples which carry 
overflow traffic, provided the grading capacity is not substantially in- 
fluenced by the network that precedes the grading. 

I. INTRODUCTION 

The Equivalent Random method 1 - is used for engineering many 
of the telephone overflow-networks in the Bell System. However, the 
method is not directly applicable to the analysis of graded-multiple 
trunk-groups t which carry overflow traffic; e.g., gradings used as 
alternate routes in step-by-step switching systems having common 
control and alternate-routing capability. Apparently, Lotze 5 -" is the 
only author with results for estimating load-loss relations for gradings 



t The reader should have some knowledge of graded multiples and the methods 
associated with the engineering of telephone overflow-networks. Some familiarity 
with the step-by-step switching system would also be helpful. Those not ac- 
quainted with these concepts may find it worthwhile to consult Ref. 1 for a dis- 
cussion of telephone overflow-networks. An introduction to graded multiples is 
given in Ref. 3. Reference 4 contains a description of the pertinent aspects of the 
step-by-step system. 
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which carry overflow traffic. Unfortunately, his method requires data 
which cannot be obtained in a step-by-step system at reasonable cost. 

We extend the Equivalent Random method to cover such applica- 
tions. The key to the extension is a technique for taking correlation 
into account when combining dependent streams of overflow traffic. 
(In principle, the technique is applicable wherever a stream of over- 
flow traffic is divided, submitted to independent trunk groups, and 
then recombined.) In Section II, we derive the appropriate covariance 
function. In order to describe how the covariance function is used to 
extend the Equivalent Random method, 1, we begin with an example 
of an application of the extended method for the analysis of a step- 
by-step graded multiple. 

Figure 1 represents schematically a step-by-step grading which 
might be used as a final route. Each horizontal bar denotes one trunk 
(server). The traffic offered to the grading* is an overflow stream 
from a subordinate network. The traffic is represented by the mean 
a and variance v of the number of simultaneous calls that would be 
in progress if this traffic were carried on a full-access group without 
blocking. The diagram is designed to indicate that an arriving call is 
first directed at random to one of the four first-choice subgroups; 
i.e., an arrival is directed to the ith subgroup with probability p { . 
After reaching a particular subgroup, the call hunts vertically upward 
for an idle trunk. If all three trunks in the subgroup are busy, the 
call overflows into the second major level of the grading. The call 
then seizes the lowest idle trunk in the second-level subgroup. If both 
trunks in the second level are busy, the call overflows to the third 
major level containing five trunks (usually called finals). If all five 
finals are busy, the call leaves the system and does not return; i.e., the 
call is blocked and cleared. The symbol a„ and v Q denote respectively 
the mean and variance of the overflow from the grading. 

In a step-by-step system, the arrangement of line-finders and 
selectors through which calls reach a grading causes an inherent load- 
balancing 4 over the first-choice subgroups ; there is positive correlation 
between the numbers of occupied trunks in the individual subgroups. 
Attempting to introduce the correlation into our model, we assume 



t The Equivalent Random method is known to be adequate for estimating 
load-service relation for overflow-networks having Poisson input. 7 - 8 

* For the present study, we assume that all offered loads are constant. Hence, 
we are estimating single-hour capacities of gradings. Utilization of our results 
for normal engineering involving average busy-hour loads, would require periph- 
eral operations to adjust the load-service relationship to reflect the effects of low, 
medium, or high day-to-day variations in the load. 9 
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Fig. 1— Schematic representation of a step-by-step graded multiple containing 
21 trunks. 

the configuration given in Fig. 2. The four arrows above the full- 
access group denote (correlated) overflow streams caused by the 
corresponding input streams. The intensity of the ith input stream 
is at = p,a. 

At this stage, we have modeled an arrival process having mean a 
and variance v. The individual substreams to the first-choice sub- 
groups of the grading are certainly correlated. How well the correla- 
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Fig. 2 — A model for correlated input streams. 
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tion approximates the correlation that exists in an actual step-by-step 
system is a point which must be tested. In Section III, we show that 
the approximation works quite well when a large group of selectors 
is connected to the grading. In any event, we can view the above 
figure as a first approximation for the step-by-step problem. More- 
over, the general features of the above configuration are not restricted 
to step-by-stcp systems. 

Figure 2 also illustrates the two basic features of the Equivalent 
Random method. First, the method assumes that, for engineering 
purposes, only the first two moments of an overflow process are 
required. Second, the method assumes that any overflow process having 
mean « and variance v is adequately approximated by the overflow 
from a unique "equivalent system" consisting of a full-access trunk- 
group with Poisson input. Whenever a and v are known, standard 
techniques are available 1 to obtain the equivalent system of c trunks 
and intensity a = a t + a-.. + a : < + a 4 of the Poisson input. 

The next step consists of using results obtained independently by 
Descloux 10 and Lotze" to determine the mean a, and the variance Vi 
of the overflow (due to a*) which is submitted to the ith first-choice 
subgroup of our grading. This "splitting" is determined by 10 

a, = PiOt (1) 

and 

y . 1 ( v 



a-LpjL-ij. (2) 

The covariance c u between the ith and jth "split" streams is given by 10 

Ca = PiP<( v — «)• ( 3 ) 

After splitting, our system is represented as 

Ca > v o) 

t 



t t t t 

a 1, V 1 a 2) V 2 fl 3' V 3 a 4' V 4 

Fig. 3 — Graded multiple with correlated input streams. 
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where it is understood that the parameters «, and v ( of the individual 
substreams satisfy equations (1) through (3). Using the Equivalent 
Random method, we now approximate the individual substreams as 
overflows from full-access trunk groups. Moreover, to determine the 
total overflow which goes to the second major level of the grading, 
we view the entire system in the following manner: 
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Fig. 4 — First cycle of application of the extended equivalent random method. 

Figure 4 indicates that an overflow of mean «, and variance v t re- 
sults from a Poisson stream of intensity a 1( being offered to a full- 
access group of Cx\ trunks, i — 1, •■•, 4. Furthermore, the ith sub- 
stream is offered to the three trunks in the ith first-choice subgroup 
of the grading, and causes an overflow of mean an and variance v n 
to be submitted to the second major level of our grading. 

One can see the main reason for looking at the grading in the manner 
described above: the parameters a,, and v u are easily computed by 1 

a,, = o»Bi 1 , i|H (On) 
and 

L (fin + 3) +on — (hi + 1J 
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where E ha (a) denotes the first Erlang loss-function (Erlang-B block- 
ing probability) . 

To complete the first cycle of computation, we need to obtain the 
mean a and variance v of the total overflow submitted to the second 
major level of our grading. The mean a is given by a = £<., a » • 
Unfortunately, the variance v is more difficult to obtain, since the 
individual substreams are correlated. 

To obtain v we need the covariance, cov (i, j), between the ith and 
j'th overflow streams which are offered to the second level. A reasonable 
method for computing cov (», j) has not been available in the past. 
Our extension of the Equivalent Random method consists of an algo- 
rithm for computing cov (t, j) in a fairly efficient fashion for many 
configurations of interest. A derivation of the algorithm is given in 
Section II. 

After cov (i, j) has been determined, v is obtained from 

4 4 

v = 2>k + £ cov(i, j). (4) 



Having a and v, we reduce the system configuration to that shown 
in Fig. 5. 

C« , v o) 



t 

(a,v) 

Fig. 5— Starting point for the second cycle of the extended equivalent random 
method. 

The proportion of the traffic (a, v) offered to the first subgroup of the 
second level is pi = (a u + an) /a, and p 2 = (a 13 + a 14 )/a is the pro- 
portion offered to the second subgroup. Consequently, one cycle of 
computation is completed. 

Repetition of the logic described above will yield estimates of the 
overflow mean a and variance v . Moreover, the cyclic nature of the 
procedure allows the logic to be programmed on a digital computer 
so that load-service tables and other relevant information can be 
generated in a straightforward manner. 
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II. MATHEMATICAL MODEL 

A service system S is composed of a collection 971 of c first-choice 
servers, two groups 9lj , 9l 2 containing ch , d 2 second-choice servers 
respectively, and two last-choice groups £ x , £ 2 each containing an 
infinite number of servers (see Fig. 6). The arrivals into the system 
are generated by v independent groups of customers G y , • • • , G, . The 
arrivals from group G, occur according to a Poisson process with in- 
tensity a, . All service times throughout the system are independent 
and have a negative-exponential distribution with unit mean. 

A customer, arriving to find an idle first-choice server, selects an 
idle server from 9TC, and service commences immediately. If an arrival 
from group { , % = 1, 2, occurs when all c of the first-choice servers 
are busy, but at least one of the d t servers from the second-choice 
group 91, is idle, a server is selected from 31, , and service commences 
at once. If a customer from group 6', arrives to find all the servers 
busy in both 9TZ and 91, , then he is served by one of the servers from 
the group £, . If fc ^ 1 and k ^ 2, requests for service from group 
G k , which occur when all c servers in 911 are busy, are dismissed and 
do not return. 

We assume that the system is in statistical equilibrium and define 
M, Ni , L { to be the number of busy servers in 911, 91, and <£, respec- 
tively (at a random instant of time) for i = 1, 2. Define the state of 
the system to be (M, N t , JV 2 , L x , L 2 ) with joint probability density 
function /(m, n v , n 2 , \ x , l 2 ) = 

P{M = m, Ni = ?h , N 2 = n 2 , L, = U , L 2 = U. 
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Fig. 6 — System configuration. 
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Setting a = a, + a 2 + ■ • « + a, , it follows that / must satisfy relations 

of the following form: 

For ^ m ^ c - 1, ^ n, ^ d h and k ^ 0, 

(a + m + rii + n* + k + Z 2 )/(m, n x , n 2 , Z, , k) 

= aj(m — 1 , ?ii ,n 2 ,li , h) 
+ (m + l)/(m + 1, w, , n 2 , l x , l 2 ) 
+ («i + l)/(™> Wi + 1, n 3 , h , k) 
+ (n 2 + l)/(m, Tii , n a + 1, h , k) 
+ (Zi + l)f(.m,ni ,n 2 , U + 1, Z„) 
+ (k + l)/(w,n, ,w 2 , U , k + 1). (5) 
Similar relations hold on the boundary of the state space { (m, n, , n-> , 
h , h) : ^ m ^ c, ^ n, ^ d«, U § 0}. We define / to be zero at all 
points not in the state space. 

The preceding infinite set of equations is quite difficult to solve. 
However, when it suffices to know the various moments of the random 
variables L x and L 2 , the problem can be simplified by introducing a 
two-dimensional binomial-moment generating-f unction. This function 
is defined by 

B(m, iiy ,n 2 ; x l , x 2 ) 

= E E Km, n 1>n2t h , l 2 )(l + x 1 )Hl + x 2 ) l ° (6) 

1,=0 f a =0 

for - 1 ^ «( £i 0, S m ^ c, and ^ th %. <h . Assuming that the 
binomial moments 



00 00 



k 

Bi I .i,(m,n ] ,n 2 ) — 2 E 



/(m, n, , n 2 , fci , fc 2 ) (7) 



exist, it follows that* 

00 00 

B(m.n x ,n 2 ;x l , x 2 ) = £ £ B«».i.(w»i **i .wO*^'' 

Of course, the binomial moments B lllla (m, n^ , n 2 ) are the entities 
of interest since 



t Various manipulations of these double series will be carried out in the sequel. 
The mathematical justification for the validity of the manipulations can be 
obtained from Ref. 11, Sections 5.3 through 5.5. 
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In particular, B<i i0) = tf(Li), £ (0 ,i) = #(£2), and B (I>1) = EiLyLo) 
so that cov(Li , L 2 ) = B (lflJ - B (1 , n) B ro ,i) • 

Relations for the binomial moments are obtained by multiplying 
both sides of Equation (5) (and the boundary equations which were 
not given) by (1 + r,)''(l + z 2 )'' and summing on /, and / 2 . Equating 
like powers of x, and x 2 yields the following finite system of equations: 

UO^m^c-1 and g /»,- £ d t , 
(a + m + n, + n 2 + /, + k)B,,.,,(m, n, , n 2 ) 

= aB tl . lt {m - l,n, , w 2 ) + (w + l)Bf,.i,(m + l,"i , n 2 ) 
+ (n, + DBL.i.Cw.n, + l,n 2 ) 

+ (n 2 + OB,,.,.^,*!, ,n 2 + 1). (9) 

ForO ^ vti ^ d, - 1, 

(a, + a 2 + c + ft, + n 2 + /, + / 2 )/?,,.,,(c, ?i, , n 2 ) 
= aBi,,,,(c — 1, n, , >j,.) 
+ a\Bi,.i~.(c, rii - 1, n 2 ) + o 2 Bt,,t,(c, n, , n 2 - 1) 
+ (», + l)B,,.,,(c, «, + l,n 2 ) 

+ (n 2 + lJB,,.,.^,*, ,n 2 + 1). (10) 

Whenever < w, ^ d, — 1, 
(a, + c + n, + d, + /, + W5i..i.(c, », , d 2 ) 

= aB,, tlt (c - l,n, , (/,.) + a,B,,,,,(c,n, - 1, d 2 ) 
+ (hB,,,,,^,^ , d 2 - 1) 

+ (n, + l)B lltll (c f n, + l,d 2 ) + a^,,.!.-^.!!, ,d 2 ). (11) 

A similar result holds for ^ n- 2 ^ d 2 - 1. At the extreme boundary 
point (c, di , d 2 ), 

(c + d 1 + d 2 + /, + WB^.j.Cc.d, ,d 2 ) 

= aBi,. t ,(c - 1, d, , d 2 ) + aiB,,,,,^, d, - 1, d 2 ) 
+ "^..Jm/, , d 2 - 1) 
+ a.B,, _,.,,(<:, d, , d 2 ) + OaB,..,.-^, d, , d 2 ). (12) 
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Since B 0| o {m, Wi , n-,) = P{M = m, N x = n x , N-, = w 8 }, it follows 
that 

c d, d, 

EE E B Q , Q (m,n v ,n 2 ) = 1. (13) 

m = n,-0 n,-0 

We set B tl ,iXm, »i > »») = ° for an y P oint & > * 2 > w > n » » n ») not in 
the set 

{(/, ,U,m,n l) n 2 ): Z, fe 0, ^ m £ c, £ n, £ d<}. 

A very useful relation can be obtained by summing all of equations 
(9) through (12) to obtain 

d, 

(Z, + Z»)B(i t ,i a) =a. E Bi,-i.i.(e, d, , n 2 ) 

n a =0 

+ a 2 t5« 1 .«,-i(c.ni ,d«). ( 14 ) 

ni -0 

Consequently, B ( i,i, = SfLiLs} can be obtained from 

-'■ 
qi(m,n 2 ) = £ B,.o(m,n, , n 2 ) (15) 



and 



That is, 



q 2 (m,n l ) = £ B 0A {m,n, ,n 2 ). (16) 



2B (ll i, = Oi? a (c, c/i) + a 2 9i(c, d a ). ( 17 ) 

Relations for gj and g 2 are obtained directly from equations (9) 
through (12) (see Appendix B). However, the relations require 
B 0f0 (c, m , d 2 ) and B , (c, d, , n 2 ) for ^ «, ^ d, . (See Refs. 12 and 
13 for a related problem.) Setting U = h = in equations (9) through 
(13) yields, with equation (14), a system of (c + 1) (di + 1) (4 + 1) 
independent linear equations for B„,o which in principle can be solved 
numerically. Unfortunately, for the step-by-step applications de- 
scribed in the introduction, c can be quite large (20 or more) although 
di and d. normally do not exceed 5. Consequently, systems of 500 
and more equations would not be uncommon. Since a solution might 
be required for several sets of parameters in any particular network, 
a direct numerical solution is not attractive. 

In Appendix A, we obtain a closed-form expression for B„, (m, ih , 
rh) in terms of the (di + 1) (da + 1) constants 
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^ w, ^ d, . (18) 



Furthermore, these constants satisfy (di + 1) (d* + 1) independent 
linear relations, so that a numerical solution is quite practical. In fact, 
the maximum size of the system of equations requiring solution for 
step-by-step systems is reduced from more than 500 to 36. 

A closed-form expression for q x (c, d 2 ) and q 3 (c, di) is derived in 
Appendix B. These results combine into the following computational 
algorithm for cov(L] , L 2 ) : First of all, 

0(0, 0) = E 1 , c (a) (19) 

(the Erlang-B blocking probability for the first-choice group), and 
for ^ rc, ^ d, , r), + n 2 > 0, 

fai + n 2 K 1+n ,/3(w 1 ,n 2 ) 

= a^w, - 1, w 2 ) + a 2 /3(n! , n 2 - 1) 



— a. 



In, - 1J 



/3(d, , n 2 ) - a 2 



d 2 

Ln 2 — 1 



Kn lf d 2 ). (20) 



By definition, /3(n, , - 1) = /?(- 1, n*) = for ^ ih ^ d, . The 

numbers v„ are intimately related to various aspects of overflow sys- 
tems 1 and satisfy the following recurrence relation: 
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for n > 1. 
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From Appendix A, 

E{Li} =» ihfi(di , 0) and E{L 2 ) = a 2 /3(0, d 2 ), (25) 

so that 

2 cov (L, , L 2 ) = a 1?2 (c, di) + a 8 gi(c, d 2 ) - 2a 1 a 2 0(di , 0)0(0, d a ). (26) 

Equations (19) through (26) constitute an algorithm for the com- 
putation of cov(Li , L 2 ). 
Whenever 

a x = a 2 and di = d 2 
the symmetry of the problem (see Fig. 1 and equation (18)) implies 
that 

/3(w! , n 2 ) = fi(n 2 , »i). (27) 

The symmetry required by (27) would prevail for most step-by-step 
graded multiples. In such cases (27) can be used to reduce the number 
of equations to M> (d, + 1) (dj + 2). Consequently, the dimensions 
of the systems of equations needed for an analysis of most step-by-step 
graded multiples would not exceed 21. Such systems can be solved 
very efficiently by numerical matrix inversion. 

III. NUMERICAL RESULTS 

In order to establish a base for comparison, we used a simulation 
to obtain load-service relations for a 25-trunk and a 45-trunk step- 
by-step graded multiple. We obtained results for several values of 
variance-to-mean ratio z = v/a and several selector configurations. 
We found that the extended Equivalent Random method furnished 
adequate estimates of blocking probability in each case where the 
maximum number of selectors was used. However, as the number of 
selectors was reduced for a particular grading, the inherent load- 
balancing 4 caused actual grading capacity to be higher than indicated 
by the extended Equivalent Random method. Consequently, we con- 
clude that the extended Equivalent Random method provides adequate 
estimates of load-service relations for graded multiples which carry 
overflow traffic, provided that the network through which calls reach 
the grading does riot significantly influence grading capacity. 

In view of the effort required to obtain the covariance cov (Li , Lj), 
one naturally questions the necessity of accounting for the correlation 
which occurs in our problem. In fact, on several occasions, we en- 
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countered the following question: What sort of results would be 
obtained if both mean and variance were split by proportion (i.e., 
aj = p [a and V\ = piV) when splitting is required, and cov(L ; , Lj) 
were assumed to be zero whenever a variance recombination is re- 
quired? 

In order to consider the preceding question, as well as to obtain a 
better understanding of the behavior of cov(L,-, Lj), three computer 
programs were written to generate load-loss relations for step-by-step 
graded multiples. These programs were based respectively on the 
following assumptions: 

(i) The input traffic is completely random, i.e., Poisson. + 
(ii) The input traffic is nonrandom. Mean and variance are split 
by proportion when required and cov (L, , L,) is assumed to 
be zero.* 
(in) The input traffic is nonrandom, and the gradings are analyzed 
by using the extended Equivalent Random method. 

Throughout the study, the offered traffic was assumed to be balanced 
over the subgroups of any grading under consideration (i.e., p t = Vi 
for all i, j) . 

For comparison, we used each of the three assumptions to compute 
load-service relations for a 25-trunk and a 45-trunk graded multiple. 
The results are displayed in Fig. 7. For these examples, the variance- 
to-mean ratio of the nonrandom offered traffic was held constant 
at 2.25. 

The different results arising from assumptions (ii) and (Hi) were 
surprising. It was originally felt by some that assumption (ii) would 
cause over-trunking, but not by the amounts indicated. For example, 
notice that the 45-trunk grading yields a B.01 blocking probability 
for an offered load of 330 ccs (9.17 erlangs) with a variance-to-mean 
ratio of 2.25 when the correlation is neglected as outlined in assump- 
tion (ii). However, Fig. 7 indicates that the same traffic can actually 
be handled at B.01 on the 25-trunk grading when the correlation is 
taken into account. Hence, assumption (ii) leads to at least 80 percent 
overtrunking for the example. An examination of other portions of 
the curves yields similar results. Consequently assumption (ii) must 
be discarded. 

Lower bounds to the load-loss relations for nonrandom traffic 



* Assumption (i) is known to cause an underprovision of trunks. 1 

* Assumption (ii) was considered by many to be the most natural approach 
to improving on assumption (i). 
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Fig. 7— Load-loss relations for the 25-trunk and the 45-trunk graded multiples. 
Variance-to-mean ratio of the offered load is z = 2.25. 



result from assumption (i) as illustrated in Fig. 7. For these examples, 
undertrunking by 18 to 25 percent results from approximation (i). 
Since the disparity will increase for larger variance-to-mean ratios, 
assumption (i) does not seem applicable either. 

Two other approximations were also tried but the results were very 
poor. The first used the correct splitting equations (1) through (3) 
but assumed cov(L ; . , Lj) = 0. As a result, some of the load-service 
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curves actually intersected each other. The second also used the 
correct splitting formulas but assumed the correlation coefficient. 

,, T v _ cov (L, , L.) 
P(L, , L t ) - [yar (L) ^ (L)]} 

to be constant at recombination points and equal to the correlation 
coefficient for the split (nonrandom) offered traffic. The predicted 
blocking resulting from the last approximation actually decreased as 
the intensity of the offered traffic increased. 

IV. CONCLUSIONS 

We have presented a technique for taking correlation into account 
when combining certain dependent streams of overflow traffic. The 
result was used to define an extension of the Equivalent Random 
method; an engineering approximation for estimating the capacities 
of overflow networks. 

The extended Equivalent Random method yields good estimates 
of load-service relations for graded multiples which carry overflow 
traffic provided the network through which calls reach the grading 
does not significantly affect grading capacity. Consequently, for ap- 
plication to step-by-stcp gradings, the extended method is restricted 
to gradings which are connected to large groups of selectors. We are 
currently investigating techniques which would allow us to consider 
systems which do influence grading capacity. 
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APPENDIX A 

System State Probabilities 

In this section, we obtain the solution for the system of equations 
(9) through (13). We use a generalization of a technique employed 
by Kosten. 14 
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For notational simplicity, let p(m, n x , n 2 ) = B 0l0 (m, n t , n 2 ), and 
let pi (m, Tii , n 2 ) denote any function which belongs to the family a 
of functions satisfying equation (9) for ^ m < oo and ^ n { ^ d, . 
Define 



Pit, ? ll92 )=EEE VM> *i , n 2 )t m qVqV. (28) 

m-0 11,-0 n,-0 

It follows from equation (9) that 



( 1 - ( )f + ( 1 -«-»| + (1 -^l = a(1 - i)P - (29) 

This linear first-order partial differential equation can be solved by 
(Lagrange's) method of characteristics (see Ref. 15, Chap. 2). The 
characteristic equations arc 

dt dq, dg 2 dP 



1 - t 1-g, 1 - tf 2 a(l - t)P 

with solutions 

\ =L ^7 = fc< and k = e~ al P 
1 t 

where k and k, are arbitrary constants. Hence, the solution to equation 

(29) is given by 



Pit, <fc , q 2 ) = ^(Y^f . TT? ) 



(30) 



where tfUi , z 2 ) denotes any analytic function of the arguments 
z x , z-> . From (28), it follows that the Taylor series for H must be 
finite, and so 

Pit, q, , q,) -f±t -(». ■ A^TiT^tT- (31) 

Following the notation of Riordan (Ref. 1, p. 89), let {a k {m): 
m = 0, 1, •••} be the sequence with generating function 
(1 — t)~ k exp(at), that is, 

The variables o»(f?») satisfy the following recurrence relations. 1 - 14 
For m s= and k ^ 0, 

ma k im) = a<T k im — 1) + fcer fc+1 (wi — 1), (33) 



CORRELATED TRAFFIC-STREAMS 

k<r k+1 (m) = (k + m — a)a k (m) + aa k ^{m), 



and 



2 a k (m) = (r k + l (j). 

m = 

It is convenient to define 

(r_,(m) = <x k (-l) = 0. 
Hence, from Equations (31) and (32), 
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(34) 

(35) 
(36) 



P(t, Qi , Q2) 



00 d i rfa 



= EEZ «fa> ,n 2 )a ni+ni (w)r(i - ?1 ) n, (i - «.r (37) 



m = ni=0 n fl =0 

00 <! i (fa 

= EEE 

m = n,=0 n,-0 



(-ir +n ' i: i: 

ii-n, fc,-n, 



ft, 

in, J 



k 

n 2 ) 



a(fc, , fc 2 ) 



•0* I+ *.(ra) 



r^'ga'. 



(38) 



Comparing equations (29) and (38), one can see that the functions 
Pi in Ct are of the form 



d ' d ' ilc 

Pi(m,n, ,n 2 ) = (-l) n,+ni £ 2 ' 



d, d, 

£ I 

A 1 = n i A- a = n a 



' l *2 
?*2 



a(&i , fc 2 )(rA 1+Aa (w), (39) 



and are determined up to the (dj + 1) (d 2 + 1) constants {a(/h , n 2 )\ 
^ n, ^ di}. There are (tf, + 1) (d 2 + 1) independent linear equations 
among equations (10) through (13), and so it follows that there is 
exactly one function, p*, in CI which satisfies equations (9) through 
(13). The appropriate restriction of p* must be p. The remainder of 
this section is devoted to a derivation of the relations which the con- 
stants {«(«, , n 2 )\ must satisfy in order to obtain the solution. 

An equivalent but less complex set of boundary conditions is ob- 
tained by putting m = c in (9) and subtracting equations (2) through 
(13) respectively. Hence, if ^ n< ^ d, — 1, 

(a — a, — a 2 )p(c, n, , n 2 ) + a,p(c, n, — 1, n 2 ) + a 2 p(c, n, ,n 2 — 1) 

- (c+ l)p(c+ l )Wl ,n 2 ), (40) 
for 0^ ?ii ^di - 1, 



2032 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1971 

(a — a r )p(c, n x , d 2 ) + a{p{c, n, — 1, d 2 ) + a 2 p(c, n x , d 2 — 1) 

= (c+ l)p(c + 1,», ,<*»), (41) 
and a similar relation holds for ^ n^ i - 1. Finally, 

ap(c, rf, , d 2 ) + «ip(c, d, - 1, d 2 ) + a 2 p(c, a", , d, — 1) 

= (c+ l)p(c + l,d, ,d 2 ). (42) 
Now, define 

P m (qx , g 2 ) = 2 Zp(m. n i. w 2 )g"'92', (43) 

ni=0 n.-0 

G...,(g2) = E p(™> n » » %)$"> C 44 ) 

n,=0 

and 

d, 

H m ,nMi) = E p(i»»«i ,w 2 )g?'. (45) 

n,-0 

It follows from equations (40) through (42) that 
[a - a,(l - ft) - a 2 (l - g 2 )]Pc(?i , fc) + ai(l - gO^Gr.d.Cft) 

+ 0,(1 - q 2 )qi'H eid Xq i ) = (c + l)P e+ ,(g, , ?,). (46) 
Equations (39) and (43) through (45) imply that 

,i, d, 
Pm(qx , q 2 ) = E E «fai ,w 2 )(T ni+na (m)(l - g,) n, (l - q 2 ) n ', (47) 

n,=0 n„=0 

<?..*.(&) = (-!)"' E «(<*i ,n,V„, +I ..(c)(l - q 2 r, (48) 

and 

ff..*.(fli) = ("I)" 1 E «(». i rf 2 )<r ni+da (c)(l - 9,)"'. (49) 

Using the identity g, = 1 — (1 — <?,) and the relations (47) through 
(49) in (46) obtains 

d, d, 

a E E afoi ,w 2 )(7 ni+na (c)(l - gi) n, (l - 92)"' 

n,-0 n,=0 



di + 1 d 

E I 

n,-l n,=0 



- O, E E «(»i ~ 1, tt a )<r n , +ni -,(c)(l " ?i) n '(l - 92)"' 
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- a 2 E E <*(»i »»a - l)(r ni+n ,_,(c)(l - g,) Bl (l - g 2 ) n * 

n, =0 n,-l 

-(-Ira."!: E (-!)"■ 



In, - 1J 

■a(rf, ,n 2 )<r,, I+n ,(c)(l - g,) n, (l - 9 2 )" a 

d 2 



- (-l)'V E E (-1)"' 



n,=0 n a = l 



n, — 1 



•a(w, , d 2 )<r ni+rfl (c)(l - g,) n, (l - q 2 )"' 
= (c + 1) E E «(«i ,fk)*lH4..(« + 1)(1 - 9,) n '(l - ffa)". 

m=0 n,=0 

Equating the coefficients of (1 — gi) n, (l — ? 2 )'" yields 
[aa ni+n ,(c) — (c + l)o- ni+n ,(c + l)]a(n, , n 2 ) 

= [aia(n { — l,n 2 ) + f^n, , n 2 — 1) + c^afa, , n 2 — l]<r ni+n ,_,.(c) 

a(d, ,n 2 )o- dl+I10 (c) 
a(w, , d 2 )<r n , +lIl (c). 



+ (-l)"> + '"a l 


f rf > 1 
7*i — 1. 


+ (-D n ' +d 'a 2 


tf 2 | 
n 2 — 1 



(*>) 



Using (33), we see that 

avm+nM - (c + l)<r» 1+ „,(c + 1) = -(n, + n 2 )<7„ 1+n , + 1 (c). (51) 
It is worthwhile to define 

0(»i , n 2 ) = (-1)"' + "'a(«i , >! 2 )o-„ 1 + „ 3 (c), and n = w, + n 2 . (52) 
Now, substitute (51) and (52) into (50) to obtain 
„ *»+i(c) 



r„(c) 



|9(ni , w 2 ) = a,/3(w, - l,n 2 ) + a 2 r5(w, , n 2 — 1) 
-J * W,n 2 ) 



— Clo 



d s 



/3(n, , rf 2 ) 



(53) 



n 2 — 1. 
for ^ rii ^ d,- and ra > 0. (Both sides vanish for n x = n 2 = 0.) 
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One additional relation is required. Using (39) , 

d, d % d, d, 

2 £ P(w,w! ,n 2 ) = 2 £ «( w i » WaK l+Bl (w) 



£ I 

n,=0 n,=0 



(-1)' 



L I 

n,=0 n,=0 

(n 

L" fcJ 
= a(0, 0)(r o (m). 
Consequently, noting relation (36) , one obtains 

d, d. 



n 2 

.'''2. 



(-!)*" 



1-0 n,=0 n a =0 



Thus, 



and 



i=EEE p(™. n i » n *) = a (°> °) S ff o(w) 

m = 

= or(0, 0)<r,(c). 
«(0, 0) = 



<r 1 (c) 



mO) = ^=EUay, 



(54) 



i.e., j8( 0f o) is the Erlang-B blocking probability (also known as the 
first Erlang loss-function) for the first-choice trunk group. 

Equations (53) and (54) completely determine (3(ni , n 2 ) for 
^ rii ^ di . Using (52) and (39 ) , we see that the state probabilities 
are given by 



p(m, ri! , n 2 ) 
Using the relation 



(55) 



m 




\m 
n. 




[fcl 




k — n 

m — n, 



it is straightforward to show that if 

e d, rfi 

A (m , „,.„,) = Z) S £ 

i=m t,=n, ti-nt 



I 

[mi 






p(», hi , fc 2 ) 
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for ^ m t== c and ^= n, ^= <7, , then an inverse relation is given by 
■p(m,n x ,n 2 ) 



d, d, 

E I 



= £ L E (- 1)'—(- 1)* — (- 1)*-— 



i 








n 2 . 



(56) 



M 

c . 








# 2 
n 2 . 



■ A(,-.A-, ,*.) 

for ^ »(. < c and ^ n, g d, . 

Consequently, comparing (55) and (56) we see that 

j8(wi i w 2 ) = A (ein ,, B ,) 



= tf 

Equations (14) and (55), imply that 

E{L X ) = Ox/3(d, ,0), 
and 

#{L 2 j = a 2 /3(0, rf 2 ). 
For the computation of /?(//, , n 3 ) using (53), the ratio 

_ <r»+i(c) 
"" «r B (c) 

is required for n ^ 1. A recursive relation for the ratio is obtained 
from (32) via 



(57) 

(58) 
(59) 



r„ + i(c) 



r„(c) 



1 + 



a g,-,(c) . 
n <7 n (c) ' 



that is, 



?i v„.,/ ?i — 



The first-order (nonlinear) algorithm is initiated with 

1 <Jn(c) 



v <7,(c) 



= E U M; 



(60) 



(61) 



i.e., the Erlang-B blocking probability for the first-choice trunk group. 
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APPENDIX B • 

A Conditional Mean 

In this appendix, a formula is obtained for the computation of f 

q (m, n) = £ B, , (m, n, , n) 

..,=0 

= E{Li | M = m, N 2 = n)P\M = m,N 2 =n\. (62) 

From equations (9) through (12), it follows that for ^ m ^ c - 1 
and ^ n g d 2 , 

(o -|- m + n + l)q(m, n) = aq(m - 1, n) 

+ (m + l)q(m + 1, w) + (n + l)g(m, n + 1) (63) 
and f or ^ n ^ tf 2 - 1, 
(a a + C + ?l + l)q(c, n) = aq(c - 1, n) + a 2 g(c, n - 1) 

+ (n + l)g(c, n + 1) + a.BoAc, d, , n). (64) 
Also, 

(c + d 2 + l)?(c, d 2 ) = ag(c - 1, d 2 ) + a 2 g(c, d 2 - 1) 

+ a y B , Q {c, di , d 2 ). (65) 

Using the methods and results presented in Appendix A, we can 
show that 

q(m, n) = £ (-1) bH W ^4^ for ^ m ^ e 



«r i+ i(c) 



where <o_i = 0, and 
(n + l)j>„+ia> n = a 2 co n -i — a 2 



and ^ n ^ d 2 , (66) 
«d, + cti/3(di , n) 



d a 

[n - lj 

for ^n ^ d 2 . (67) 

In particular, gCc, d->) = «>rf 2 , and can be obtained from (67) in the 
following closed form: 



t Throughout this appendix, the subscript 1 on <?i has been omitted for nota- 
tional simplicity. 
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> *) =- - 'trr . fztr r (68) 



d a 



nfr 



A similar expression is valid for g^(c, di). 
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